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Some of the recent developments in the theory of random surfaces and simplicial quantum gravity is reviewed. 
For 2d quantum gravity this includes the failure of Regge calculus, our improved understanding of the c > 1 
regime, some surprises for q-state Potts models with g > 4, attempts to use renormalization group techniques, 
new critical behavior of random surface models with extrinsic curvature and improved algorithms. For simplicial 
quantum gravity in higher dimensions it includes a discussion of the exponential entropy bound needed for the 
models to be well defined, the question of "computational ergodicity" and the question of how to extract continuum 
behavior from the lattice simulations. 



1. INTRODUCTION 

The first attempt to discretize gravity in a co- 
ordinate free way dates back to the seminal work 
of Regge from 1961 [1]. The purpose of Regge 
was purely classical. He wanted to approximate 
a manifold with a given metric by a piecewise 
linear manifold. Regge made the important ob- 
servation that parallel transport and the integral 
of curvature have a natural definition on piece- 
wise linear manifolds. For a given abstract tri- 
angulation an assignment of length to the links, 
consistent with the triangle inequalities, will de- 
fine uniquely the concept of curvature and paral- 
lel transport in a way which coincides with an em- 
bedding of the triangulation as a piecewise linear 
manifold in some ambient space R" . The assign- 
ment of integrated curvature is as follows: Let T 
be the triangulation, let {rii} denote the set of i- 
simplexes in T. A given (d — 2)-simplex nd-2 will 
belong to a number o{nd-2) of rf-simplexes, nd{k), 
k = 1,2, ...,o{nd-2), where o{nd-2) is called the 
order of nd-2- For one of the nd{k)'s we have that 
nd-2 belongs to precisely two (rf— l)-subsimplexes 
and if we denote the angle between these two sub- 
simplexes 6^, the deficit angle e{nd-2) is defined 
as 



The angle 6^ is defined in the ambient space R" 
but can be expressed entirely in terms of the 
length of the links of the abstract triangulation. 
In the same way the volume of Vi{ni) of an i- 
simplex in the triangulation can be expressed en- 
tirely in terms of the length of the links. The 
fundamental definitions of Regge, which is in ac- 
cordance with the concept of parallel transport 
on the piecewise linear manifold is given by: 
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The idea of Regge was to consider the link 
length as the classical dynamical variables. For 
a given manifold one should choose a suitable tri- 
angulation and by varying the length of the links 
find the extremum of the Einstein action. The 
corresponding piecewise linear manifold should 
then qualify as a approximation to the underly- 
ing smooth manifold which is a solution the Ein- 
stein equations. This makes perfect sense at a 
classical level since the lengths of the links are 
invariants and for a sufficiently fine triangulation 
one expects a good approximation to the smooth 
manifold. 

Inspired by lattice gauge theories the quantum 
analogue to Regge calculus was suggested. The 
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simplest idea is to use the classical link variables 
and integrate over these variables. However, one 
encounter a complicated counting problem. In 
quantum gravity we want to integrate over equiv- 
alence classes of metrics. For a given abstract tri- 
angulation there can be many length assignments 
which correspond the same equivalence class of 
metrics as is clearly seen by considering trian- 
gulations of the plane with the usual Euclidean 
metric. For a given triangulation of the plane 
there is usually plenty of room for moving the 
vertices around thereby changing the length of 
the links without changing the Euclidean metric 
of the plane. We conclude that the replacement 

involves as highly non-trivial jacobian ■J{li). In 
eq. (4) the functional integration on the Ihs is 
over equivalence classes of triangulations, sym- 
bolized by the division by the "volume" of the 
diffeomorphism group. The integration on the 
rhs of eq. (4) is over all link length compatible 
with the triangle inequalities. We can view the 
triangle inequality as incorporated in the Jaco- 
bian ■J{li). An important aspect of the "quantum 
version" of the Regge calculus is the absence of a 
cut-off. As formulated above one has discretized 
gravity but not introduced a cut-off. In general 
one will have to introduce a smallest distance or 
higher derivative terms in the discretized action 
in order to make the integration over link length 
well defined. 

A different approach was advocated by Wein- 
garten, both for strings and for gravity. The idea 
was to consider manifolds embedded in hypercu- 
bic lattices. In this case we have a natural cut 
off, the lattice spacing a, and we have no problem 
with overcounting since each manifold is explic- 
itly given on lattice as a "physical" , coordinate in- 
dependent manifold. The metric assignment for 
such a piecewise linear surface can be made by 
the Regge calculus. Since all building blocks are 
identical it will only depend on the number of d- 
and [d — 2)-dimensional building blocks of the d 
dimensional manifold on the _D-dimensional hy- 



percubic lattice. In this way we get: 
Seh[G,\] = j d'i^-^ j d'i^R^ 

S[kd-2,kd] = kdNd - kd-'jNd-'j (5) 
and the partition function will be given by 
Z[kd-2,kd]= J2 e-^['=--'=^]. (6) 

MGZO 

From this point of view the hypercubic manifolds 
form a grid in the space of equivalence classes 
of metrics. The question is to which extent it 
is uniform and dense in the limit of infinite Nd- 
While a dependence on the dimension D of the 
hypercubic lattice is natural in the case of string 
theories, it is not desirable in the case of gravity. 
Weingarten suggested a formulation independent 
of the underlying hypercubic lattice, but the more 
streamlined formulation in terms of abstract so- 
called dynamical triangulations were first given 
in two dimensions in [2-4] and in higher dimen- 
sions in [13,6,7]. In this formulation the summa- 
tion over hypercubic manifolds is replaced by the 
summation over abstract rf-dimensional triangu- 
lations where the length of the links is taken to 
be a. In this way we still have a cut off which is 
"reparametrization invariant" and by Regge cal- 
culus we can still assign a metric and the action 
(5) to the given abstract triangulation. Again the 
action is only a function of Nd-2 and Nd since 
all fc-simplexes in the triangulation are identical. 
Each triangulation is in this way a representative 
of a whole equivalence class of metrics. There will 
be no problem with overcounting and by summing 
over all abstract triangulations we again lay out a 
grid in the space of equivalence classes of metrics. 
Eq. (6) will be replaced by 

Z[kd-2, kd] = J2 6-'=^^^+'=-=^-= (7) 
TeT 

Eq. (7) describes a discretization of pure gravity. 
T denotes an abstract triangulation in a given 
class of triangulations T. In addition one can cou- 
ple matter to gravity in a natural way [8]. Since 
the link length is taken to be one it is easy to show 
that a Gaussian scalar field couples as follows: 

S[^,g]= I d'^^g-'da^db^^ 

JM 
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(8) 



(U)GT 



In this formula i and j denote vertices and (ij) 
links in the triangulation T. This formula can 
be derived by assuming that the field in the in- 
terior is given by the linear extrapolation of the 
field value at the vertices. Space-time in the 
interior of the rf-simplexes is assumed fiat and is 
dictated by the length assignment a to all links. 
Alternatively one can place the discretized field 
values in the center of the rf-simplexes and again 
arrive at an expression like (8), only is the in- 
dex i now a rf-simplex and the summation is over 
neighboring rf-simplexes. Spin systems and even 
gauge fields can also be coupled to gravity in this 
formalism. 

It is possible to divide the study of simplicial 
quantum gravity in two classes: the study of two- 
dimensional quantum gravity and higher dimen- 
sional gravity. The theory of two-dimensional 
gravity is well defined if we restrict the topol- 
ogy of the manifold. It also includes the theory 
of bosonic strings in since such strings can 
be viewed as two-dimensional gravity coupled to 
D scalar fields. It is also possible to move in 
the direction of the theory of membranes since 
one can introduce extrinsic curvature terms in 
such models. For the two-dimensional theories 
the situation is very nice from the point of view 
of computer simulations: We have systems where 
all our methods: Monte Carlo simulations, finite 
size scaling etc. seem to work well. This is illus- 
trated in fig. 1 and fig. 2. The results presented 
are from a Monte Carlo simulation of the Ising 
model coupled to two-dimensional gravity. The 
simplest finite size scaling tells us that 

Mat ~ xiv ~ AT^/'^''^ 



(9) 



where N is the volume (i.e. in this case the num- 
ber of triangles in the triangulation), Mjv the 
magnetization and xn the susceptibility at the 
pseudo critical point corresponding to N . dn de- 
notes the Hausdorff dimension of two-dimensional 
gravity coupled to a c = 1/2 conformal field the- 
ory. From eq. (9) we get 
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Figure 1. The magnetic exponent as a function 
of volume as defined by eq. (10). 



y 1.4 



.2- 



0.8 



0.6 



o 0.9196 
□ = 0.7733 



(b) 



O CD 



CD m 



10 



100 



1000 



10000 



Figure 2. The magnetic susceptibility exponent 
as a function of volume as defined by eq. (10). 



A nice convergence to the known theoretical value 
is seen. For details I refer to [9]. 

In addition there are regions, most notably 
matter fields with central charge c > 1 where 
the theory is not yet understood from an analyt- 
ical point of view. The numerical results might 
serve as an important inspiration for analytical 
approaches which so far have been successful for 
c < 1. If we discuss the higher dimensional grav- 
ity theories we move into unchartered territory 
and the study of these theories is at present at 
the explorative stage. 

In the rest of this review I will concentrate on 
the developments during the last year. For the 
two-dimensional theories this includes 
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1. The failure of two-dimensional quantum 
Regge gravity. 

2. Understanding large c. 

3. Surprises for large q-state Potts models. 

4. Attempts to formulate renormalization 
group methods. 

5. New algorithms 

For three- and four-dimensional gravity it in- 
cludes 

1. The exponential bound on the entropy. 

2. The question of "computational ergodic- 
ity". 

3. How to extract physics from the Monte 
Carlo simulations. 

2. TWO-DIMENSIONAL GRAVITY 
AND RANDOM SURFACES 

2.1. The failure of quantum Regge gravity 

While Regge 's approach to gravity makes per- 
fect sense in the classical context in which it was 
invented, it has already been mentioned that use 
of the formalism in a quantum context involves 
a non-trivial Jacobian J{li) (see eq. (4)) which 
is presently unknown. Until now it has been as- 
sumed, appealing to universality, that the choice 
of this measure is not of uttermost importance. 
A scale invariant measure 

n' 

has been used in most simulations. During the 
last year two large scale simulations have tested 
whether Regge gravity in 2d agrees with contin- 
uum calculations and simplicial quantum gravity. 
In both cases the result is so far negative.. 

The first work to be discussed is by Holm and 
Janke ([10] and the contribution to Lattice94). 
They have looked at an Ising model coupled to 
Regge gravity. They used toroidal topology and a 
regular lattice with dynamical length assignment 
of the links. The spin was located at the vertices. 
The measure was given by eq. (11) and the lattice 



sizes in the Monte Carlo simulations ranged from 
6^ to 512^. In addition to the Einstein-Hilbert 
action a term 

a J v^i?2 (^2) 

was added to the action for a = 0, 0.1 and 0.001. 
The result of a very careful finite size scaling anal- 
ysis of the critical exponents of the Ising model is 
summarized in table 1. a is the exponent for the 
specific heat, /3 and 7 was defined above, while S 
is the exponent defined by M ~ at the crit- 
ical point. Finally i> is the exponent which gov- 
erns the divergence of the correlation length as 
one approaches the critical point. It is clear from 
table 1 that the agreement with the standard On- 
sager exponents for the Ising model on a regular 
lattice is perfect, while there is disagreement with 
the known exponents for dynamical triangulated 
surfaces. Since these latter exponents agrees with 
Liouville theory one is let to conclude that Regge 
gravity coupled to Ising models does not describe 
c = 1/2-conformal matter coupled to gravity. It 
is an interesting question whether it is a funda- 
mental fiaw in the approach or it can be repaired 
by a different choice of measure (11). 

The second attempt to test the viability of 
Regge gravity in the quantum regime is by Bock 
and Vink (see [11] and contribution to Lattice94). 
They have performed extensive Monte Carlo sim- 
ulations of the two-dimensional Regge version of 
quantum gravity. They have included in addi- 
tion to the Einstein-Hilbert term a i?^ as in (12). 
In this way it is possible by careful scaling argu- 
ments to derive a relation between the coupling 
constant a in (12) and 7^ , the string susceptibility 
exponent. One definition of 7^ is as an entropy 
exponent in the partition function. The partition 
function for a given area A is defined by 

In this formula /j^c is a non-universal cut-off while 
the subleading correction A^''~^ defines the uni- 
versal critical exponent. In pure gravity it can by 
shown, either by the use of Liouville theory or by 
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Table 1 



Critical exponents in the Regge-Ising model 
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the DTS-model, i.e. simplicial quantum gravity, 
that 

7. = 2 - \x, (14) 

where x is the Euler characteristic of the surface. 
While 7s is readily identified in simplicial quan- 
tum gravity one has to use scaling arguments in- 
volving the i?^ term in the Regge formulation. 
The problem with the measure is here highlighted 
as it turns out that 7^ seems to depend on choice 
of measure: If (11) is modified to 

i h 

there will be a dependence on r\. The interpre- 
tation of this is not clear. From universality and 
the principle of ultra-locality one would not ex- 
pect a change like (15) to affect any continuum 
quantity. To be conservative Bock and Vink did 
their simulations with ?y = in order to be in 
agreement with earlier measurements of Hamber 
where agreement with (14) have been reported. 
However the improved simulations of Bock and 
Vink for the sphere and the bi-torus (x = 2 and 
—2) showed no relation between the Regge theory 
and the correct results given by eq. (14). Again 
the conclusion is that we know presently no mea- 
sure (i//(/) of Regge gravity which can produce the 
correct continuum results. 

2.2. Matter with c>l coupled to 2d gravity 

In the context of statistical mechanics it is easy 
to couple matter with a central charge c > 1 to 
two-dimensional gravity: Take the statistical sys- 
tem, e.g. a multiple Ising system. It is defined on 
a regular lattice where it has a critical point and 
a second order transition and an associated diver- 
gent correlation length. The central charge and 
the associated critical exponents are defined at 



the critical point. The operational steps needed 
in order to couple the system to 2d quantum grav- 
ity is as follows: Define the statistical system on 
an arbitrary random lattice which represents a 
triangulation (or more generally a polygon net) 
on a surface of fixed topology^. Define a new 
statistical system by taking the annealed average 
over the selected class of random triangulations 
of the surface. The annealed system will still 
have a critical point. The order of the transition 
can change (if it change it is usually to a higher 
order transition, which is intuitively reasonable) 
and at the new critical point there will be new 
critical exponents. From the continuum point of 
view we attribute this change to the interaction 
with gravity. However, also purely "geometrical" 
exponents like 7^ can change when we perform 
the summation over matter configurations and af- 
terwards the summation over triangulations with 
the effective action in part coming from the usual 
gravity action and in part coming from the inte- 
gration or summation over matter fields. We view 
this as the back reaction of the quantized matter 
on gravity. 

For c < 1 many of the annealed lattice models 
have been solved explicitly and, as already men- 
tioned above, the critical exponents agree with 
the continuum calculations using Liouville the- 
ory. However, the continuum approach breaks 
down for c > 1. This is illustrated by the so- 
called KPZ formula for 7^ : 

^ c - 1 - ./{c- l)(c- 25) , , 

74X = 2) = (16) 

It should be noticed that the formula only makes 
sense for c < 1 and that 7s (c) < for these val- 
ues of c. It has been a longstanding question what 

^Of course the definition should be reasonable, i.e. the 
part of the action which involves nearest neighbor inter- 
action should be in accordance with the geometry repre- 
sented by the random lattice. 
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happens for c > 1. While it is not clear how to 
proceed in the continuum approach, the theory 
is clearly well defined and unambiguously defined 
in the statistical mechanics approach. During the 
last year there has been an increased understand- 
ing of the interaction between gravity and matter 
with c > 1. In particular the limit c — cxd is 
now well understood by mean field calculations 
[12-14]. 

Let me briefiy describe the simplest derivation 
of the c — cxD results: First I appeal to the Monte 
Carlo results which has been obtained over the 
last couple of years: They showed that for the 
annealed average of many Ising models over ran- 
dom triangulations the following is true: 

(1) : There is a critical (inverse) temperature j3c 

and a phase transition between a magne- 
tized phase (for large [3) and a phase where 
the total magnetization is zero (for small 
/?). 

(2) : For [3 > j3c the exponent 7^ jumps to —1/2, 

the value for pure gravity. For increasing 
c, i.e. increasing number of Ising models 
coupled to gravity, the system is increasing 
fast completely magnetized just above j3c 
which increases logarithmically with c. 

(3) : At j3c one has 7s (/3c) > for sufficiently 

large c. 

(4) : In a region below j3c there seems to be a 

very strong interaction between that matter 
system and the spin system. The random 
surfaces seem strongly deformed compared 
to typical surfaces present in computer sim- 
ulations without many spin systems. Again 
the effects are enhanced with increasing c. 

(5) : For /3 — the exponent 7^ is equal —1/2, 

the value of pure gravity. 

Points (1) and (2) motivate the mean field ap- 
proximation for large c: Since magnetization is 
the same all the way from [3 = 00 io j3c the spin 
excitations will just be a selfconsistent iteration 
of the minimal energy excitation, i.e. just the 
first term in the usual weak coupling expansion. 
The major difference between the situation on a 



dynamical lattice and a fixed regular lattice is 
that the minimal energy spin excitation on a reg- 
ular lattice involves the fiip of a single spin since 
this creates the smallest boundary. Not so on a 
dynamical lattice. Arbitrary large spin clusters 
can be separated by small boundaries. This is 
shown in fig. 3 and after having realized this it 




Figure 3. Graphical representation of the mean 
field equation (17). 



is also possible to understand intuitively the in- 
teraction between spin and geometry: The spm 
excitations which will dominate for large [3 and 
large c will he the ones which deform the geom- 
etry into baby universes of alternating spins. If 
we consider the one-loop function where the loop 
just consists of two links we have the selfconsis- 
tent equation shown in fig. 3: 

G(//,/3) = ^e-''^-(l + e-2/'G(//,/3))^- 

T 

= ^e-'^^- =Go(//), (17) 

T 

where jj = /j — |log(l -|- e'"^!^ G(ji, j3)) since we 
have that the number of links Lt = f for a 
triangulation. This implies that we can solve for 
the one-loop function G{iJ,,l3) in terms of the a 
know function: Goifi)'- 

G{^,(]) = Go{-^) (18) 
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// = //+-log(l + e-2/'Go(//)) 



(19) 



The eqs. (18) and (19) are easily solved [13]. The 
solution for 7^ {j3) is shown in fig. 4. 
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Figure 4. The solution to eqs. (18) and (19) for 

Is- 



The conclusion is that we indeed can obtain 
7s > 0. In fact one gets 7^ = 1/3 at the transi- 
tion point l3c and 7^ = 1/2 below /3s . It cannot 
be correct for /3 — according to point 4) above, 
but it is seen that the mean field approximation 
explains all other observations. In particular it 
is tempting to associate the region just below /3c 
with a phase where the interaction between mat- 
ter and geometry is so strong that the typical sur- 
face consists of a number of small baby universes 
which are individually magnetized but where the 
total magnetization is zero, i.e. a branched poly- 
mer phase. The value 7s = 1/2 is precisely the 
generic value for branched polymers and a surface 
build of many small baby universes would be such 
a branched polymer surface. For small /3 the ap- 
proximation breaks down since there can be many 
other spin excitations than the selfconstistent it- 
erations of the minimal excitation discussed here. 

Is it an accident that the first value of 7s > 
which is different from the branched polymer 



value is 1/3 instead of 1/2? The answer is no! 
There exists by now a theorem which states that 
for multiple spin systems coupled to gravity one 
must have 



7s > ^ 7s 



1 



r 



n > 0, 



(20) 



for a spin system. The proof is an elaboration 
of the eqs. (18) and (19) [15] which leads to the 
following representation of any 7s > in a system 
with local spin interactions: 



Is 



Is 



Is 



1' 



7. < 0. 



(21) 



In fact the interpretation of (20) and (21) is that 
the typical surface with gg = 1/n -|- 1 will consist 
of baby universes with 7s = —l/n (i.e. c < 1) 
glued together along small loops. The generic 
value will be n = 2 which corresponds to c = in 
which case the interpretation is that the individ- 
ual baby universes are completely magnetized, i.e. 
corresponds to 7s = —1/2, and the new critical 
behavior arises at the transition where the spin of 
individual baby universes are aligned. It should 
be mentioned that the mean field equations above 
can be solved in an external magnetic field and 
one can verify that the transition is a third order 
transition. 

An interesting and yet unanswered question 
concerns the transition from the branched poly- 
mer phase below /3c where 7s = 1/2 and back to 
the phase where 7s = —1/2, which is valid for 
small 7s . 

Another interesting question is how universal 
this behavior is and if there are examples of non- 
trivial 7s >. By non-trivial I now mean a 7s dif- 
ferent from 1/2 and 1/3. In fact it is not too diffi- 
cult to find a meanfield model which is a general- 
ization the model above and which has 7s = 1/4 
[16]. It is however very interesting that a com- 
pletely different class of random surface theories 
seems to have the same critical behavior. This is 
random surface theories defined on dynamical lat- 
tices with the usuall Gaussian action and in addi- 
tion with an extrinsic curvature term. These have 
some interest also membrane physics since they 
can serve as toy models of biological membranes. 
It has also been conjectured that they are related 
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to superstrings. It is known that they have a 
phase transion for D = 3, D being the dimension 
of target space, for a finite value of the extrin- 
sic curvature coupling constant [17]. Below the 
transition point we have the bosonic string which 
has a non-scaling string tension [18]. At the crit- 
ical point it seems that the mass and the string 
tension scales and the random surface model is 
a candidate for an interesting continuum string 
theory [19,20]. This has been extensively cov- 
ered at the former lattice conferences. Using baby 
universe counting it is possible to measure 7^ at 
the transition point[21]: it seems as if gs = 1/4, 
i.e. we have an absolutely non-trivial realization 
of eq. (20). The results of the measurement of 
7s at the critical extrinsic curvature coupling is 
7s = 0.27 ± 0.02. This is illustrated in fig. 5. It 
is interesting that the simulations seem to agree 
with old simulations of on a hyper-cubic lattice 
where one can also add and extrinsic curvature 
term [22]. 

0.5 1 , , , , , , 1 

0.45 - 
0.4 - 
0.35 - 




0.15 
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0.05 - 

I ' ' ' ' ' ' 

20 40 60 80 100 120 a_min 



Figure 5. The value of 7^ for random surfaces 
with extrinsic curvature coupling at the critical 
extrinsic curvature coupling constant, amin de- 
notes the minimal area of the baby universes used 
in the determination of 7^ . I refer to [?] for de- 
tails about the determination of 7^ from the the 
baby universe distributions. 



2.3. q state Potts models, q>4 

At first it seems that these models should not 
be interesting since the models on a regular lat- 
tice have a first order transition, and from this 
point of view they do not fit into the framework 
outlined above. However, after coupling to grav- 
ity one can prove that large q = large c [12] and 
in addition it is known that the order of the tran- 
sition is often decreases when the annealed aver- 
age over triangulations is performed. For a single 
Ising model it is known that the transition is de- 
creases from a second order transition on a regu- 
lar lattice to a third order transition on dynam- 
ical lattices. Similarly, since large q corresponds 
to large c we know from the arguments presented 
in the last section that q ^ co also corresponds 
to a third order transition. What is the situation 
for finite g > 4? 

The results of extensive numerical simulations 
performed for the annealed average of g = 10 and 
q = 200 models is shown in the second and third 
row of table 2 (for details I refer to [9] and the 
contribution of Thorleifsson at Lattice94). The 
first row shows the known analytical results for 
g = 4, which corresponds to a c = 1 conformal 
field theory, while the last row shows the mean 
field results for c = q = co. 

The exponents suggest strongly that g = 10 be- 
longs to the same universality class as g = 4 while 
g = 200 belongs the same universality class as 
q = CO. This implies that 7^ = 1/3 at the transi- 
tion for g = 200 and according the the mean field 
picture there should be a region below /3c where 
7s = 1/2. The corresponding curve is shown in 
fig. 6 and gives partial confirmation of this picture 
despite the fact that it is known to be quite diffi- 
cult to extract a reliable 7s from the spin models 
coupled to gravity. An even more convincing con- 
firmation comes from a measurement of the spin 
clusters. Let 

W = ^ (.2, 

denote the ratio between the number of boundary 
links hi between spin clusters and the number Nc 
of clusters. W = 1 only if there is a total magneti- 
zation (one cluster) or if there is a tree structure 
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Table 2 

Critical exponents for the g-state Potts models. 
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Figure 6. 7s for q = 200. 

of spin clusters with only one-link boundaries^. 
The measurements of W shows that for q = 200 
W is close to 1 in a region below j3c. In the same 
[3 region Nc is large. This shows that in a re- 
gion below j3c there are many small spin clusters. 
However, the clusters with 6; = 1 is in 1-1 corre- 
spondence with the underlying geometry of baby 
universes since we can only have a spin cluster 
with boundary length one if there is a baby uni- 
verse separated from the rest of the triangulation 
by a 1-loop. If we increase [3 and pass j3c W will 
stay very close to one, but Nc will jump from be- 
ing a fraction of the total number of links to one, 
corresponding to the situation with total magne- 
tization. For q = 200 this is in perfect agree- 
ment with the mean field picture. For g = 10 
there is a clear deviation: for [3 < j3c W > 1 
and it drops to = 1 only after j3c. This is even 
more pronounced at g = 4 and on regular lattices. 
W for various Potts model on dynamical triangu- 
lated surfaces are shown in fig. 7. To summarize: 
There seems to be a critical gc > 10 such that for 

^In the simulations we have allowed degenerate triangula- 
tions with 2-loops and 1-loops. 



-<^q = 4 
-^q = 10 
-^q = 50 
-^q = 200 - 



-i -0.5 6 'o'.5 

T 

Figure 7. Plots of W for various Potts models in 
the neighborhood of their critical points. 



q < qc the q-state Potts model on dynamical tri- 
angulated lattices belongs to the same universal- 
ity class as g = 4, i.e. they describe a c = 1 theory 
coupled to gravity even if the models on a regular 
lattice cannot be assigned a central charge. For 
g > gc the models at the critical point is identi- 
cal to the meanfield model described in the last 
section. 

It is interesting that the same picture has been 
suggested for models with c > 1 in [23]. It was 
suggested that there was a region 1 < c < cq 
where the models had the same back-reaction on 
gravity as the c = 1 model. 

2.4. Renormalization group approaches 

Since finite size scaling seems to work well in 
two-dimensional gravity coupled to matter, it is 
natural to try to use another standard technique 
like the renormalization group. This technique 
might help us to define the scaling limit in sit- 
uations where we have no analytical results to 
guide our choice of scaling. If one wants to ap- 
ply the renormalization group equations there are 
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two basic steps: 

(1) define large geometrical cells. 

(2) define blocls spins or bloclsed matter fields 
on the larger geometrical cells. 

Usually (1) is considered trivial and all the inge- 
nious efforts enter into step (2). When we dis- 
cuss dynamical triangulations (2) becomes triv- 
ial compared to (1). Somehow we have to find 
in a consistent way a method of forming blocks 
which preserves the underlying fractal structure 
of space-time. An old suggestion which goes all 
the way back to the study of hypercubic random 
surface models [24] is conceptually clear: Con- 
sider dynamical triangulations in two dimensions. 
The fundamental lattice length is the link length. 
It appears in the fractal structure precisely by 
the appearance of baby-universes where the neck- 
length is minimal, i.e. one if we allow one-loops 
in the triangulations. There will be a distribu- 
tion of such baby universes as a function of their 
area. Compare this distribution to one where the 
class of triangulations does not allow one-loops, 
but does allow two-loops. Effectively this corre- 
sponds to a doubling of the cut-off length. In the- 
ory one can continue this blocking and the contin- 
uum fractal structure should be defined from the 
distribution of baby universes which represents a 
fixed point in the blocking procedure. From a 
practical point of view this procedure is useless. 
In [25] (see also contribution to this conference) 
it was noticed that the above idea could be dras- 
tically simplified if one stayed with a fixed neck- 
length of the baby universes, e.g. the minimal 
one. The distribution of baby universes of area 
B, for a fixed area A of the total two-dimensional 
universe, is given by: 

~yl5^'-2. (23) 

However, the baby universes have a hierarchic 
structure in the case of spherical topology: One 
can cut them away, starting from the outer baby 
universes, and at each step close the boundaries. 
In this way the initial surface of area A will lead 
to a sequence of ensembles of surfaces with areas 
A < Aq, the starting area: 

S{A - Ao) ^ Pi{A\Ao) ^ P2{A\Ao) ^ • • • (24) 



The distributions Pi, P2, . . . are readily available 
from the simulations and this allows us to deter- 
mine 7s under the hypothesis that the distribu- 
tion of baby universes are given by (23). 

A completely different approach has been sug- 
gested in [26,27]. The idea is to perform a block- 
ing of a random lattice in such a way that some 
important geometric features are approximately 
preserved. Start out with a regular lattice and 
perform a blocking to a coarser lattice. Again we 
restrict ourself to 2d lattices for simplicity. The 
updating of the lattice structure is performed on 
the fine lattice. Mark the vertices of the coarser 
lattice. These vertices have an existence on the 
fine lattice. During the Monte Carlo simulation 
they will "fioat" around. How to update the 
coarse lattice, i.e. how to change neighbor as- 
signment. In [26] it was suggested to let the link- 
distance on the underlying fine lattice determine 
whether or not to fiip a link on the coarse lattice. 
This is illustrated on fig. 8: Two triangles ahc and 
acd have the link ac as common link. While the 
vertices ahcd have an existence on the finer lat- 
tice the links ah, ac etc. only exist on the coarser 
lattice. The standard move on a two-dimensional 
triangulation is a link fiip: ac hd. In order 
to decide if one should perform this fiip the link 
distance between a and c and between h and d is 
calculated on the underlying fine lattice. If it is 
smaller between h and d the link is fiipped from 
ac to hd. The blocking of matter fields which are 
assigned to the vertices is straight forward. 

In [27] it was shown that the prescription works 
well in pure gravity. If one adds formally irrel- 
evant operators like i?^ to the Lagrangian and 
perform the blocking there is a clear fiow to pure 
gravity. 

It is important to make the above methods into 
genuine quantitative methods which allow us to 
extract critical exponents and correlation length 
in a reliable way. 

2.5. Algorithmic improvements 

In any Monte Carlo simulations the efficiency 
is intimately linked to the importance sampling. 
Whenever one encounter the phenomena of crit- 
ical slowing down it is important to choose a 
sampling of configurations which is in accordance 
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Figure 8. The flip of a link in the coarse grained 
lattice. The link distance between a and c in the 
underlying fine lattice is 3, while it is only 2 be- 
tween h and d. 



with the critical fiuctuations. If the nature of 
the critical fiuctuations is understood it might be 
possible to diminish the autocorrelation time be- 
tween successive configurations by orders of mag- 
nitude. The best known example is the cluster 
algorithms for spin systems. These can in fact 
be used without any modifications on dynamical 
triangulations. For large central charge c of the 
matter fields coupled to 2d quantum gravity the 
fractal structure becomes very pronounced as de- 
scribed above. The fiip which is used to move 
around in the class of triangulations with fixed 
area becomes quite inefficient. This makes the 
computer simulation of large systems impossible. 
Recently it has been shown how the autocorre- 
lation time for the dynamical triangulations can 
be drastically improved by using so-called baby 
universe surgery [28]. In the case of pure two- 
dimensional gravity one can cut away any baby 
universe outgrowth which is connected to the rest 
of the universe by a minimal bottle-neck and glue 
it back at any other position at the surface. This 
operation will not change the action at all, but if 
the baby universe is large it will be a global move 
which can only be imitated by many local moves. 
In addition it is perfect in accordance with the 
pronounced fractal structure of the surfaces for 
large c. In this way it has a lot in common with 
the cluster algorithm for spin systems. Using this 
algorithm it will be possible in the future to sim- 
ulate systems of order 10® triangles. This would 



be very time consuming relying entirely on local 
moves like the fiip. 

There is clearly a lot of room for algorithmic 
improvement of the Monte Carlo simulations of 
dynamical triangulations and the baby universe 
surgery algorithm is presumably just the first of 
many improved algorithms. 

3. QUANTUM GRAVITY FOR d>2 

3.1. Definition of the model 

Could it be that the formalism of simplicial 
gravity, which works so well for 2d gravity can 
serve as a definition for quantum gravity in c? = 4? 
We are clearly moving into unchartered territory 
since there exists presently no continuum theory 
of quantum gravity in c? = 4. The Euclidean path 
integral is not particularly well defined: First the 
Einstein-Hilbert action is unbounded from below: 
For any given metric assignment to a manifold, 
gab{x), a conformal transformation 

gab{x) Q{x)^gab{x) (25) 

can make R{x) arbitrary positive if only the 
derivatives of il{x) are sufficiently large. In addi- 
tion we have to address the question of summa- 
tion over topologies in the path integral 

For c? < 3 there is no discussion of the meaning of 
Top. The concept of homeomorphism is in one- 
to-one correspondence with the concept of diffeo- 
morphism. However, for c? > 3 this correspon- 
dence breaks down: there are topological struc- 
tures which allow infinitely many smooth struc- 
tures and there are topological structures which 
do not allow any smooth structures. 

If we decide that the X^lxop means the sum- 
mation over smooth structures it is possible to 
replace "smooth" by "piecewise linear" or "com- 
binatorial inequivalent" . For d < 7 there there 
is a one-to-one correspondence between "piece- 
wise linear" and smooth structures. In this way 
simplicial gravity as described in the introduction 
offers a unique possibility to write eq. (26) in a 
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united way: 

Z{k2,ki) = ^e-'=^^^(^)+'==^=(^). (27) 

T 

where the summation is over all combinatorial in- 
equivalent abstract triangulations. 

3.2. The entropy bound for fixed topology 

Although the definition (27) is tantalizing from 
the point of view that it unites topology and met- 
ric structures in a natural way, it is unfortunately 
formal. This is true even in c? = 2. The sum is di- 
vergent. A limiting procedure should be taken in 
order for sum to make sense. In two dimensions 
the so-called double scaling Imitt is an attempt to 
perform the summation. 

In c? = 4 it is not yet known how to perform the 
summation. However, it is generally believed that 
eq. (27) makes sense if one fixes the topology and 
the first numerical simulations seemed to support 
this idea [29]. A necessity is that the number of 
abstract triangulations of a given combinatorial 
manifold is exponentially bounded as a function 
of #4, the number of four-simplexes. This is seen 
in the following way: If the topology is fixed it 
is easy to show that for any triangulation T with 
N4{T) large we have 

l#4(T)^#2(T)^y#4(T). (28) 

This implies that Z{k2, k^) defined by (27) exists 
for a fixed topology if and only if there exists a 
^4 such that 

Mop(A^4)= J2 <e'^'^'- (29) 

Tjv4Gtop 

For a given ^2 the lowest possible ^4(^2) is the 
critical value of the bare cosmological constant 
^4. 

Until now the only topology which has been 
used in the computer simulations is that of 5"* and 
the bound (29) has recently been questioned [30]. 
However, successive simulations clearly favor an 
exponentially bound on the number of triangula- 
tions [31,32]. Very recently an analytical proof of 
the exponential bound has been published [33]. 



3.3. Computational ergodicity 

It is known that combinatorial four-manifolds 
are not algorithmically classifiable. This makes it 
questionable that one will ever be able to simu- 
lation the complete sum (27) numerically even if 
we manage to define it in one way or another. It 
is also known that there exists four-dimensional 
manifolds which are not algorithmically recogniz- 
able in the class of all four-dimensional combi- 
natorial manifolds. Does this imply anything for 
Monte carlo simulations on such manifolds? This 
question was answered affirmatively two years ago 
by Ben-Av [34]. 

First a few definitions: Two triangulations are 
called combinatorial equivalent if there exists a 
triangulation which is a common subdivision of 
the two. The Monte Carlo simulations use a local 
update of the triangulations, which only requires 
a finite, fixed number of operations for each step. 
The updates are called moves. In addition the 
finite set of moves used are ergodtc: Given two 
combinatorial equivalent triangulations Ti and T2 
one can get from one to the other in a finite num- 
ber of steps. 

In case we consider triangulations of a manifold 
Aio which is not computational recognizable one 
can prove that if 

#4(^1) < N and #4(72) < N (30) 

the number of steps m the finite algorithm to get 
from any triangulation Ti of A4o to any other 
triangulation T2 of A4o cannot be bounded by a 
recursive definable function r{N). This includes 
functions like N\ and N\^'. For a given N the 
total number of triangulations with #4(T) < N 
is clearly bounded by some {aN)\ where a is of 
order one. How can the above mentioned theo- 
rem be true? It can only be true if we on the 
way from some Ti to some other T2 is forced to 
triangulations with very large N4. In this way 
there will be barriers separating different regions 
of configuration space. In fact one can prove that 
the height of the barrier cannot be bounded either 
by a recursively definable function. 

Although in principle ergodic, the class of lo- 
cal moves mentioned above cannot in practise 
reach all triangulations of A^o- However, it is not 
known if the region of configuration space we can- 
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not reach is of measure zero in the limit N ^ co. 

How is the situation for 5"* , which is the topol- 
ogy used until now. It is not known if 5"* is algo- 
rithmically recognizable or not. Extensive com- 
puter simulations looking for any sign of barriers 
have failed (see [35] or the contribution to Lattice 
94). This can be taken as a hint that either 5"* 
is computational recognizable or a hint that the 
regions which cannot be reached are of measure 
zero. It would be most interesting to perform the 
simulations for a manifold which is known not to 
be algorithmically recognizable. 

3.4. Scaling 

The following scenario seems to be true in four- 
dimensional simplicial gravity: For large values of 
the bare gravitational coupling constant (small 
values of ^2) the typical universes are very crum- 
pled with almost no extension if distance is mea- 
sured as the shortest link length between two ver- 
tices. For small values of the gravitational cou- 
pling constant (large values of ^2) the universes 
seem very elongated, and the Hausdorff dimen- 
sion seems close to two: the universes seem to be 
branched polymers!. There is a phase transition, 
most likely of second or higher order where the 
Hausdorff dimension might be close to four!. 

The different groups which have performed 
computer simulations agree more or less on this 
scenario. The next step is to analyze in detail 
the scaling behavior in the infinite volume limit 
close to the phase transition. The first step has 
been taken in [36] where it was found that dis- 
tribution of internal geodesic distances seems to 
allow a scaling which indeed indicated a fiow to- 
wards the transition point in the following way: 
The distributions of distances for universes corre- 
sponding coupling constant away from the critical 
point could be mapped on distributions for large 
universes close to the transition point. This is the 
first genuine hint of a divergent correlation length 
when we approach the critical point. 

Clearly much remain to be done. In particu- 
larly it would be helpful with simulations of larger 
universes. Until now most simulations have used 
N4 < 32.000. Although a large number such N4S 
correspond in fact to quite small four-dimensional 
manifolds. 



4. Summary 

Two-dimensional gravity is a wonderful play- 
ground for both numerical and analytical work. 
Many of the tools we have available from sta- 
tistical mechanics seems to work well and some- 
times the models formulated on dynamical trian- 
gulations are much more easily solved analytically 
than the counter parts defined on a regular lat- 
tice. This might be a refiection of the underlying 
reparametrization invariance of the theory in the 
scaling limit. 

Higher dimensional gravity is still unchartered 
territory. Hopefully the numerical simulations 
will teach us something about the typical quan- 
tum universe and inspire analytic work. Dynam- 
ical triangulations have proven a very powerful 
tool in two dimensions and one can hope that the 
same will be true in higher dimensional gravity. 

There is a lot of room for algorithmic improve- 
ments in the simulations. The first step has been 
taken by the use of baby universe surgery. It 
would be very helpful if we could increase the 
size of the simulated universes by some order of 
magnitudes. 
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